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Abstract 

Molecular center-of-mass motion is taken into account in the theory of co- 
herent control of nuclear spin isomers of molecules. It is shown that infrared 
radiation resonant to the molecular rovibrational transition can substantially 
enrich nuclear spin isomers and speed up their conversion rate. 
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I. INTRODUCTION 



Symmetrical molecules exist in the form of nuclear spin isomers [[[]]. Most famous are 
the ortho and para isomers of H 2 . @. On the other hand, spin isomers of molecules heavier 
than hydrogen are almost unknown. This is because their investigations were not possible 
for a long time due to the lack of enrichment methods. Now such methods start to appear 
that has advanced the field significantly (see the review ||). Yet one needs more efficient 
enrichment methods in order to perform investigations of spin isomers of various molecules. 
It will also stimulate practical applications of spin isomers. 

It was proposed recently to use strong electromagnetic fields for the spin isomer enrich- 
ment. In the first version of the method an enrichment is produced by the light-induced 
change of molecular level populations MM- In the second method enrichment is due to 
combined action of the light-induced level populations change, level energy shift and light- 
induced coherences The latter allows to refer to this method as coherent control of 
spin isomers. 

Analysis of the isomer coherent control ||[7| was performed using simplified model of 
molecules interacting with a strong electromagnetic field. Molecules were assumed to be at 
rest and have nondegenerate states. The purpose of the present paper is to develop more 
general model which will account molecular center-of-mass motion. This can be essential 
because the process takes place in the gas phase at room temperature where light-molecular 
interaction is known to be affected by the molecular motion. 

II. QUALITATIVE PICTURE 

Coherent control of spin isomers has been analysed for the case of microwave excitation of 
molecular rotational transitions || and for infrared excitation of rovibrational transitions |7| . 
Although each of these two schemes has some advantages we will account molecular motion 
only for the case of rovibrational excitation because this scheme allows to achieve larger 
enrichment. Note also, that molecular motion is more important for infrared excitation 
than for microwave excitation because the latter has smaller Doppler broadening. 

Let us start with the qualitative description of the phenomenon. We assume that the 
molecules have two nuclear spin states, ortho and para. This is the case, for example, for the 
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CH3F molecules |]|| • Nuclear spin states are different by the symmetry of nuclear spin wave 
function of equivalent nuclei in the molecule. Long wavelength electromagnetic radiation 
cannot change the nuclear spin state directly. An important point is that the spin conversion 
in polyatomic molecules is governed by intramolecular mixing of close (accidentally resonant) 
ortho and para states (see Ref. Jl for more details). Interaction of an external radiation 
with these particular states can affect the spin conversion. 

The level scheme for rovibrational excitation is shown in Fig. [I]. Similar to Ref. J7| we 
assume that each of the two vibrational states (ground and excited) has one ortho-para level 
pair mixed by intramolecular perturbation. The infrared radiation excites the rovibrational 
transition m — noi the ortho molecules. By assumption, collisions cannot change spin state 
of the test molecules directly, i.e., collisional cross-section of ortho-para transitions vanish. 
Suppose that the molecule is placed initially in the ortho subspace of the ground vibrational 
state. Due to collisions the test molecule will undergo fast rotational relaxation inside the 
ortho subspace. This will proceed until the molecule jumps to the state m', which is mixed 
with the para state k! by the intramolecular perturbation V' (see Fig. [l]), or to the state 
n which is mixed with the para state k by the combined action of the external field and 
intramolecular perturbation V. Admixture of a para state implies that the next collision 
can move the molecule to another para states and thus localize it inside the para subspace. 
The spin conversion occurs. 

Vibrational relaxation in this enrichment scheme plays an important role. After ortho- 
para transition in upper vibrational state the molecule relaxes to the ground vibrational 
state. Consequently there is no back conversion through the m — k channel in excited 
vibrational state but back conversion occurs through the unperturbed m'—k' channel. Steady 
state concentrations of the spin isomers is determined by the ortho-to-para conversion rate in 
the excited vibrational state and the para-to-ortho conversion rate in the ground vibrational 
state. It is essential that conversion through the m — k channel can be substantially enhanced 
by electromagnetic radiation. 

It is clear that molecular center-of-mass motion affects the phenomenon because it makes 
the ortho-para mixing in upper vibrational state velocity selective. One obvious effect is that 
the radiation does not interact with all molecules in the state n. Molecular motion influences 
also the level energy shift and coherences produced in the molecule by external radiation. 
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One needs quantitative account of these effects in order to understand the real efficiency of 
the method. 



III. RATE EQUATIONS 

Quantitative analysis of the problem can be performed with the help of kinetic equation 
for the density matrix, p. The molecular Hamiltonian consists of four terms, 

H = H + hG + hV + KV'. (1) 

The main part, Ho, has the eigen ortho and para states shown in Fig. |l|. TiG describes 
the molecular interactions with the external radiation that will be taken in the form of 
monochromatic travelling wave, 

G = -(E d/^)cos(cj L t-kr), (2) 

where E ; ojl and k are the amplitude, frequency and wave vector of the electromagnetic 
wave, respectively; d is the operator of the molecular electric dipole moment. V and V' are 
the intramolecular perturbations that mix the ortho and para states in excited and ground 
vibrational states, respectively. These perturbations originate from the nuclear spin-spin, or 
spin-rotation interactions P-TT . 



In the representation of the eigen states of H kinetic equation reads [12| , 



dp/dt + vVp = S-i[G + V + V',p], (3) 

where S is the collision integral; v is the molecular center-of-mass velocity. We have consid- 
ered molecular motion classically in this equation and neglected spontaneous decay which is 
slow for molecular rovibrational transitions in comparison with collisional transitions. Ac- 
count of molecular motion makes the matrices p and S velocity dependent. The matrices 
G, V, and V remain to be velocity independent. 

Kinetic equation for the total concentration of molecules in one spin state, e.g., ortho 
can be obtained from Eq. (|3|), 

dp /dt = 2Re / i(p m kV km + Pm'k'V^ m ,)dv; Po = J2 P(i> v ) rfv ; 3 G ortho. (4) 

3 J 

Here j indicates the rotational states of the molecule. These sates are assumed to have no 
degeneracy. In Eq. (H) a uniform spatial distribution of molecular density has been assumed. 



Collision integral did not enter into Eq. (||) because by assumption collisions do not change 
the molecular spin state, i.e., J2j I Sjjdv = 0, if j G ortho, or j £ para. G did not enter into 
Eq. (|4]) also because the matrix elements of G off-diagonal in nuclear spin states vanish. 

The off-diagonal matrix elements p m k and p m >k' can be found using the perturbation 
theory. We assume further the perturbations V and V being small and consider zeroth- 
and first-order terms of the density matrix, 

p = p<°> + p« (5) 

Collisions in our system will be described by the model standard in the theory of light- 
molecular interaction. The off-diagonal elements of S have only decay terms, 

s 3f = - r Pii'; j f- (6) 

The decoherence rates, T, were taken equal for all off-diagonal elements of collision integral. 
Relaxation of the off-diagonal elements of the density matrix is determined in molecular 
systems mainly by the level-population decay which is thus assumed to be j-independent. 
Note that there are opposite examples, e.g., [|T3|]. 

The diagonal elements of S can be expressed through the kernel of collision integral, A, 
as follows, 

S( a ) — ^2 J A(a\cti)p(ai)dvi — p(a) ^2 J A(ai\a)dv 1 ; a = {a,j, v}. (7) 

Here a = e,g denotes the excited and ground vibrational states, respectively. The nuclear 
spin index is omitted in Eq. (^) because there are no collisional ortho-para transitions and 
the kernels A for ortho and para states are assumed to be the same. We will consider further 
the model of strong collisions with the following collision kernel, 

= v r w{])5 aai 5{\ - vi) + v t w{j)f{v)5 aai + v v w(j)f{v)5 aie 8 ag . (8) 

Here w(j) is the Boltzmann distribution of rotational state populations, 

w(j) = Z- 1 exp(-E J /kT) } (9) 

where Z is the rotational partition function for one vibrational state; Ej is the rotational 
energy; T is the gas temperature; is the Boltzmann constant. We assume the same 
function w(j) for the two vibrational states of both isomers, /(v) is the Maxwell distribution, 
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f(y) = ~^ ex P(-^ v Q = pk B T/M, (10) 

7T ' Vq V 

where M is the molecular mass. 

In Eq. (|8|) the three relaxation rates were introduced, rotational relaxation (z/ r ) without 
change of molecular velocity (velocity persistence effect) and vibrational state; translational 
relaxation (u t ) that equilibrate also rotational distribution but do not change the vibrational 
state; vibrational relaxation (u v ) which is accompanied by velocity and rotational relaxations. 
The introduction of different relaxation rates for different degrees of freedom makes the 
model of strong collisions more accurate and flexible. It allows to adjust the model to the 
particular experimental conditions. Because of its simplicity model of strong collisions is 
often used in laser physics, see, e.g., fH||r5[ . 



IV. ZEROTH-ORDER THEORY 

For the zeroth-order term of the density matrix one has the following kinetic equation, 

dp^/dt + v ■ Vp (0) = S {0) - i[G, p (0) ]. (11) 

Electromagnetic field excites only ortho molecules. Consequently, para molecules remain at 
equilibrium in the zeroth-order perturbation theory. For the population of para molecules 
one has, 

/fG^v) = (iV-pi°V(j)/(v); pf (e,j,v) = 0, (12) 

where N is the total concentration of the test molecules; pf> is the concentration of ortho 
molecules. 

Zeroth-order theory for ortho molecules is in fact a familiar phenomenon of light- 
molecular interaction. There are numerous publications on this subject including many 
textbooks. Our zeroth-order theory is close to the approach jlSj , although the relaxation 
scheme given by Eq. (H) is different from [|1^]. Eqs. an d fllTf ) allow to deduce the 

system of equations for the stationary populations of ortho molecules, 

u eP V>{e,j, v) = iy r w(j)p^(e, v) + v t w(j)f(v)$\e) + p^p(Q, v) 5 jm ; 
u gP ®{gJ,v) = v r wU)P { o\9^) + vM3)f(v)p { °\g) + 
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Pi°tov) = r2 + 2 ( ^ 2 kv)2 K 0) G?,n,v) -pW( e , m ,v)] , (13) 

where i/ e = z/ r + z/ t + v v \ v g = v r + v t \ p^(e) and p^(g) are the concentrations of ortho 
molecules in excited and ground vibrational states; p^(a, v) = pi°H a 5 J5 v ); a = 
and p(f2, v) is the probability of optical excitation. Eq. ([13|) is written in the rotating wave 
approximation. Nonzero matrix element of G is given by, 

G mn = -Ge^""*); G = E d mn /2h, (14) 

where f2 = u;^ — i^ mn is the radiation frequency detuning from the absorption line center, 
co> mn ; the bar over a symbol indicates a time-independent factor; |G| is the Rabi frequency. 

Solution of Eqs. (|H|) has no difficulty Integration of Eqs. (JL3|) over v and summation 
over j gives the system of algebraic equations for the distributions p^(a,j, v), p^(a,v), 
and p^ (a) that can be easily solved. The level populations can be expressed as the sum of 
three terms, equilibrium in j and v, equilibrium in j, and nonequilibrium in j and v, 

pW(e,j, v) = p(°){p(0)r ev «;(i)/(v) + v)[r et ™(j) + r e 5 jm ]} 

pPfo.J.v) = p?> Mi)/(v)[l-p(0)r s „] -K^v)[r gtW (i) + V in ]}. (15) 

here p(fi) = / p(f2, v)dv and, 

ro , w(n)T%f(v) 

nii ' vj ntrl + ^-kv^fi + i/K + F^Vn]- 1 J 

In this equation, Tb = Ty/1 + k is the power dependent homogeneous linewidth (Bennett 
width) and k = 2|G| 2 rir _1 is the saturation parameter. The relaxation times are expressed 
through the collisional parameters of the kernel A from Eq. (|8|), 

T ev = v~ 1 -(vy + vt)' 1 ; r et = (u v + vtY 1 - is' 1 ] r e = v~ x 
Ti = T et w(m) + r gt w(n) + T e + T g ; r 2 = T ev w(m) + T gv w(n). (17) 

The function, 

Y(Q) = ^yReW{x + iy); x = Q/kv ] y = T B /kv , (18) 
is the Foigt profile of the absorption coefficient. W(z) is the probability integral, 
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1 + — / e u du 
n Jo 



(19) 



Similar to the calculation of the diagonal matrix elements of pf\ one can find the off- 
diagonal density matrix element, p^(n\m; v) which amplitude is equal to, 

rt°>(n|m;v) = "^ 0) ^V - <(fl - kv)]. (20) 



V. FIRST ORDER THEORY 

In zeroth-order perturbation theory, one neglects perturbations V and V' . It implies that 
there are no coherences between ortho and para states, p^ k = 0; p^, k , = 0. Consequently, 
one has, 

dp /dt = 2Re J i(p^ k V km + p%V^ t )dv, (21) 

instead of Eq. (|J). Thus the spin conversion appears in the second order approximation. 
The first order correction to the density matrix, is determined by the equation, 

dp m /dt + v- Vp (1) = S (1) - i[G, - i[V + V, p(°>]. (22) 

For p^, k , one has from this equation, 

& = ^B\pf\9,k'^) -pW(g,m',v)\, (23) 

J. — r~ LUJ 

where uj' = oj m /y. The density matrix element p^ k can be obtained from equations which 
are deduced from Eq. fl22"|). 

d Pmkl m + V • V Pml + T Pmk + ^mnPifc = iV mkP ( o ] {e, TTl, v); 

dp^/dt + v • Vp2 + TpS + zG nm p« = iV mkP W(n\m; v). (24) 

Substitutions, V mk = Ve iujt , (to = to mk ); p ( ^ k = p^ ] k e iut ; p { 2 = ptl^^ 1 '^, transform 
Eqs. ( p^) to algebraic equations from which one finds p^ k . Using this result for p^ k and 
p^) k , from Eq. (p3|) one has the rate equation ( ^T|) in the form, 

'p?(9,k')-p^(g,m') - 



d Po _ 2T\V'\ 2 
~dt ~ T 2 + lu 
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o\ V \2 Re I [T + i(n + LJ- kv)]p(°)(e, m, v) + tGp^(n\m; v) ^ ^ 



(r + iu)[r + i(n + uj- kv)] + \g\ 2 

This rate equation has the same general structure as in [0. On the other hand the contri- 
butions from the level population in upper state, p^(e, m, v), and from the light-induced 
coherence, p^(n\m; v), are more complicated now due to their velocity dependencies. 
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VI. ENRICHMENT AND CONVERSION 



Similar to Ref. @], the rate equation (|25f) is presented in the form, 

dpo/dt = Ni op - p Q r, 7 = lop + l'po - In + 7n + Icoh- (26) 

We have neglected in the right-hand side of Eq. fl26|) the small difference between p^> and 
the total concentration of ortho molecules p Q . In Eq. ( PS] ) the following partial conversion 
rates were introduced. The field independent rates, 

2iW'| 2 2r|^'| 2 
< P = r2 + ^/2 w ( fc/ ); ipo = T 2 + U ,2 w ( m ')- ( 27 ) 

f ree = j' op + r y' po , is the isomer conversion rate without an external field. The rates 7' and 
7p G coincide with the analogous rates in pj. This is not surprising because intramolecular 
ortho-para state mixing is velocity independent. The field dependent contribution to the 
conversion rate, 

l'n = l'poP(ty W ( m ')( T 9V + Tgt), (28) 

is also similar to that in |7j except for the new definition of relaxation parameters. This 
rate enters to the total rate 7 in Eq. (|26| ) with minus sign because it is caused by the light- 
induced population depletion of the ground vibrational state of ortho molecules. The last 
two contributions to 7 are more complicated now. The "noncoherent" term appears due to 
the level population, p^(e,m, v), in Eq. (pop, 



[r + i(Q + uj - kv)}[p(Q)T ev w(m)f(v) + p(Q,v)(T et w(m) + r e )] 



(T + iu))[r + + uj — kv)J + \G\ 



The "coherent" term, •y co h, originates from p^(n|m;v), in Eq. fl25|) , 

^ = ! lr+ *^-* a -?* *. (30) 

J [F + iu))[r + z(i2 + uj — kvjj + \G\ Z 

General expressions for the conversion rates are rather difficult to analyse. It will be 
done elsewhere. Below we compare the results [0] with the present calculations. Solution to 
Eq. ( Hp) can be presented as, 

Po = p + (Po(0) -p )exp(-7t); p D = ^7^/7. (31) 

Here p Q is the stationary concentration of ortho molecules. Without an external radiation 
(at the instant t = 0), the equilibrium concentration of para molecules reads, 
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P M = N- Po (0) = iV 7 ; o /7/ree = N/2, (32) 

if the Boltzmann factors are assumed to be equal, w(k') = w(m'). This implies 7' = 7^ 
(see Eq. (^7])). A stationary enrichment of para molecules is defined as, 

^Jj)- 1 = 1 - 2 T' (33) 

One can see from this equation that the external field produces an enrichment of para isomers 
if 7 ^ 7 /ree . 

We assume in further analysis the following parameters, the Doppler width of the 
absorbing transition, kv = 30 MHz; uo = 100 MHz, J = 130 MHz, V^, k , = 5 kHz, 
Y/P — 2 • 10 8 s _1 /Torr and the Boltzmann factors of the states m', k! , m, and k all equal 
10 -2 . This set of parameters gives the field free conversion rate, 7/ re e = 10 -2 s _1 /Torr, which 
coincides with the conversion rate in 13 CHsF. Nuclear spin conversion in these molecules is 
governed by quantum relaxation (see the review [[|). The translational, rotational and 
vibrational relaxation rates will be taken equal, v t = 0.1Y, v T = 0.1Y, and v v = 0.01T, 
respectively. 



Close similarity between noncoherent, 7 n , and coherent, 7 co h, terms (see Eqs. ( ^9|) and 
(p0|)) suggests that their contributions to the total rate, 7, are rather difficult to separate. 
It is indeed so for a wide range of parameters. On the other hand, at some conditions it is 
possible. For example, it can be done for weak external fields, thus small G. In this case one 
has, similar to 0], two peaks in enrichment at frequencies Q ~ —uo and Q = 0, see Fig. |2|. An 
obvious difference is that these two peaks are much wider now having the Doppler width, 
kvo. The data shown in this figure correspond to V m k = 3 kHz, and T = 2 MHz. The peak 
at Q = appears because the excitation probability, p{Q), has maximum at Q = 0. This 
peak has positive contribution from the rate 7 n and negative contributions from the rates 
7^ and 7 C oh- At high G the main contribution originated from 7^ reaches the value, 

In ^ 1 ( 2 T ev w(m) + T et w(m) + T e 

7 f ree 2 \V'u J T\ + T 2 

which constitutes ~ 15% in our numerical example. The homogeneous width of the peak at 
Q = is equal to Yb- Consequently, this peak has rather large power broadening. The peak 
at Q = in isomer enrichment was predicted in ||] by considering only the level population 
effects. 
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The "off-resonant" peak at ~ —uj appears because the external field splits the ortho 
state m into two components and one sublevel crosses the para state k at this frequency of 
the external field (see Ref. for more details). At low G this peak is determined by 7^. 
When G increases the peak at Q ~ — uj grows to much bigger values than the peak at Q — 0. 
At high G the main contribution to the Q ~ —uj peak comes from the rate j n . Position 
of this peak is given by — u[l — \G\ 2 /(T 2 + uj 2 )], and the homogeneous width is equal to, 
T[l + |G| 2 /(r 2 + uj 2 )], thus this peak has much smaller homogeneous width than the peak 
at VI = if uj > r. 

The data shown in Fig. correspond to strong optical field, G = 50 MHz, and three 
values of V m k- T was taken equal 2 MHz. One can see, that strong optical field is able to 
convert almost all molecules to the para state if V m k — V m 'k'- Thus relatively weak (3 kHz) 
coupling in upper state produces macroscopic effect, viz., almost complete enrichment of spin 
isomers. In the case of much weaker coupling in upper state, enrichment is still significant. 
For example, if the perturbation in upper state, V m k = 100 Hz, one has the enrichment, 
/3 ~ 1%. Enrichment at this level can easily be measured. The enrichment peak at Q ~ —uj 
is narrow (having almost the Doppler width) and thus can be distinguished from much wider 
structures (the width ~ T B ) caused by the = peaks. 

Conversion rate in the system is equal to 7, see Eq. (|3ll) . It is convenient to characterize 
it in relative units, 

Irel = l/l free ~ 1- (34) 

Thus, jrei = without an external field. At low G, conversion rate, like enrichment, has 
two peaks in its frequency dependence. If Rabi frequency, \G\, is large and the ortho-para 
couplings in upper and low states have the same order of magnitude, conversion can be 
significantly enhanced (Fig. |3|, upper panel), although this enhancement is smaller than in 
the case of the molecules at rest 0. This enhancement appears because of the crossing 
of ortho and para states in upper vibrational state by external field. The main difference 
between the present case and Ref. is that the peaks in conversion rate are much wider 
now having the Doppler width, kvo- 
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VII. CONCLUSIONS 



We have developed the model of coherent control of nuclear spin isomers taking into 
account the molecular center-of-mass motion. This model is more realistic than the model [[/J 
in which the molecules were assumed to be at rest. It was shown that the main conclusions 
of |7j still hold. One has two peaks in the enrichment and conversion rate dependences 
on the laser radiation frequency. One can achieve a substantial enrichment and speed up 
the conversion rate. It is also valid that this enrichment scheme can be used to access 
very weak perturbations in molecules. On the other hand, molecular motion causes some 
decrease in sensitivity of the method. For example, 1% enrichment can be achieved for 
the intramolecular perturbation equal V ~ 100 Hz in the present model, in comparison 
with 1% enrichment at V ~ 30 Hz in Ref. 0. This is not a significant difference. But 
even this loss in sensitivity can probably be avoided in special excitation scheme, e.g., by 
using a multifrequency excitation instead of a monochromatic excitation. This needs further 
investigations. 

Search of proper resonances will be one of the main difficulty in practical implementation 
of the proposed enrichment scheme. An existence of the off-resonance peak in the isomer 
enrichment has particular importance. First, one does not need to search for an exact 
coincidence of the ortho and para states in upper vibrational state. The energy gap can 
be compensated by detuning the laser frequency. Second, the excitation in the off-resonant 
peak does produce even higher enrichment than in the resonant peak but does not produce 
significant population in upper vibrational state. The latter helps not to populate the 
excited vibrational state of other spin isomer and consequently suppress the spurious back 
conversion through the upper channel. 
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FIG. 1. Level scheme. Vertical solid line indicates infrared excitation. Vibrational relaxation is 
shown by the dashed vertical lines. V and V indicate the intramolecular ortho-para state mixing. 
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FIG. 2. Frequency dependence of the enrichment of para molecules, (3, at the Rabi frequencies 
\G\ = 1, 2, and 4 MHz. 
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FIG. 3. Enrichment of para molecules, (3, and conversion rate, j re i, for the three values ofV m k- 
The wide negative structures in the middle and low panels (the width equal Tb) are mainly due 
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